On the geometry of classically integrable 
two-dimensional non-linear sigma models 

N. Mohammedi * 

Laboratoire de Mathematiques et Physique Theorique ( CNRS - UMR 6083), 
Universite Frangois Rabelais de Tours, 
Faculte des Sciences et Techniques, 
Pare de Grandmont, F-37200 Tours, France. 



Abstract 

A master equation expressing the zero curvature representation of the equations of 
motion of a two-dimensional non-linear sigma models is found. The geometrical prop- 
erties of this equation are outlined. Special attention is paid to those representations 
possessing a spectral parameter. Furthermore, a closer connection between integrabil- 
ity and T-duality transformations is emphasised. Finally, new integrable non-linear 
sigma models are found and all their corresponding Lax pairs depend on a spectral 
parameter. 
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1 Introduction 



The problem of finding dynamical systems which are integrable is a fascinating subject in 
mathematics and theoretical physics. In classical mechanics integrability is understood as 
the possibility of finding as many conserved quantities as the number of degrees of freedom 
of the dynamical system. It happens, in some cases, that these conserved quantities lead 
to the exact solvability of the associated equations of motion. In field theory, however, an 
infinite number of conserved charges is required for integrability. 

The Lax formulation of integrability provides a method for constructing conserved dy- 
namical quantities. In this formulation, a two-dimensional field theory is considered to be 
classically integrable if a Lax pair (A ,Ai) can be found such that the linear syste nfl 

[d Q + Ao(\)]V = 

[<9i + Ai (A)] * = (1.1) 

yields, as its consistency condition, the equations of motion of the two-dimensional theory 
under consideration. Here the matrices Ao and A\ depend on the fields of the theory and 
possibly on some free arbitrary parameter A, known as the spectral parameter. This pa- 
rameter can be very useful in extracting conserved quantities. The fields \l/ can be either 
a column vector or a matrix of the same dimension as Ao and A\. The consistency con- 
dition (usually referred to as the zero curvature condition) of this linear system is clearly 
{d A 1 -d 1 Ao + [A ,Ai]}V = 0. 

The conserved quantities are then constructed using the so-called monodromy matrix 

T(A,r) = Pexp (- J** Ai (A, a, r) da) , (1.2) 

where P stands for the path-ordered exponential and we have chosen a to be in the interval 
[0 , 27t]. One can show that the traces of powers of the monodromy matrix, Tr [T n (A, r)], 
are independent of the time r and are in involution with respect to Poisson brackets: 
{Tr [T m (Ai, r)] , Tr [T n (A2, t)]} = 0. The proof of the first statement assumes the peri- 
odicity condition Ao (A , , r) = Ao (A , 2n , r). Expanding Tr [T n (A, r)] in powers of A 
generates an infinite set of conserved charges (see [TJ [2] for more details). 

In this paper we would like to examine the question of integrability in two-dimensional 
non-linear sigma models. This is because there are only a handful cases of such theories which 
are known to be integrable (the principal chiral model, the Wess-Zumino-Witten model and 
their various modifications [3j HJ El E]). It is therefore important to investigate whether 
other integrable models exist. Furthermore, the study of the properties of non-linear sigma 
models involves often the geometry of the target space on which these theories are defined. 
For instance, the renormalisation properties of these models constrains the geometry of the 
target space [7]. It will be shown in this paper that the requirement of integrability puts 
further constraints on the allowed target spaces. This could be of crucial importance to 
string theory as non-linear sigma models are supposed to describe the propagation of the 
massless modes of bosonic string theory [8]. In other words, the conditions for conformal 

1 Here, the two-dimensional coordinates are (t, a) with do — and d\ — In the rest of the paper, 
however, we will use the complex coordinates (z = r + ia , z — r — ia) together with d = and d — -IL. 
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invariance at the quantum level (the vanishing of the beta functions) and the requirement 
of classical integrability of non-linear sigma models might reduce the number of possibilities 
for the spaces on which one can carry out the compactification of the extra dimensions of 
string theory. 

We start this paper by giving the general framework for a zero curvature representation 
of the equations of motion of a two-dimensional non-linear sigma model. We derive a target 
space condition for this requirement and analyse its resulting geometry. In section 3, we 
provide some known solutions to this condition. Section 4 deals with the issue of introducing 
a spectral parameter in the Lax pair construction and further geometrical properties are 
analysed there. We then study, in section 5, the integrability of a non-linear sigma model 
which generalises the principal chiral sigma model. In section 6, the interplay between T- 
duality and integrability of non-linear sigma models is explored. This work is a continuation 
of an earlier investigation [5] . 

2 Zero curvature representation of non-linear sigma 
models 

A two-dimensional non-linear sigma model is an interacting theory for some scalar fields 
tp l (z , z) as described by the action 



S 



J dzdzQij ((p)d<p i d(p i . (2.1) 



The metric and the anti-symmetric tensor fields of this theory are defined as 

9ij = 2 V^V Qji) i ^ij — — [Qij Qji) ■ (2-2) 

We will assume that the metric is invertible and its inverse is denoted g % K Indices 
are raised and lowered using this metric. We will also define, respectively, the Christoffel 
symbols, the torsion and the generalised connection as follows 

r| = ^g kl (digij + djgii - %y) 
= ^(dibij + djbu + dibjt) 

= Y%-H% . (2.3) 

The equations of motion of this theory can be written as 

S l = Bd<p l + tydip'dcpi = . (2.4) 

Let us now construct a linear system whose consistency conditions are equivalent to these 
equations of motion (a zero curvature representation). We take, as an ansatz, this linear 
system to have the following form 



* = 



<9 + /3, (<^)<V * = , (2.5) 
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where CKj and ft are two matrices depending on the fields ip % . This form of the Lax pair 
is dictated by the fact that the equations of motions of the non-linear sigma model do not 
contain terms involving d 2 or B 2 . 

The compatibility condition of the linear system takes then the form 

T = (ft - at) ddip 1 + {dtfj - d j0li + [a, , &]) dtp'dtpj = . (2.6) 

The non-linear sigma model enjoys a zero curvature representation of its equations of motion 
if this compatibility condition can be written as 

jr = = (2.7) 

for some matrices fa {ip). In order for this last relation to yield S 1 = as the only non trivial 
possibility, the matrices fa have to be linearly independent and their number must be equal 
to the dimension of the target space of the non-linear sigma model. 

The compatibility condition of the linear system yields the equations of motion of the 
two-dimensional non-linear sigma model, that is equation (12. 7p holds, provided that the 
matrices (</?), ft {ip) and fa {ip) satisfyH 

ft OC{ fa 

difr - d 3 ai + [on , ft] = fa . (2.8) 

The first equation gives simply ft in terms of a, and fa 

ft = «i + fa . (2.9) 

The second equation of the above set can then be written as 

F {j = - (Vifa - ft£ ;fx k ) , (2.10) 

where we have introduced, for later use, the field strength Fy and the gauge covariant 
derivative corresponding to the matrices a» 

WiX = diX+[a l} X] , (2.11) 

where X denotes any matrix valued quantity. 

Equation ( 12 .101) is at the centre of the ability to represent the equations of motion of 
a non-linear sigma model as a zero curvature condition of a linear system. The unknowns 
of the problem are the two sets of matrices on and fa and the generalised connection Q^. 

Each triplet (ati , fa, fi^J satisfying (12.101) . yields a non-linear sigma model with a zero 
curvature representation (provided that one can extract gij and bij from the knowledge of 
f&). However, equation (12.10p does not guarantee that the matrices on and fa will depend 
on a spectral parameter (which plays an important role in the construction of the conserved 
quantities of two-dimensional integrable theories). Let us now explore some properties of 
this central equation. 

2 These equations are found by comparing the terms involving ddip 1 and dip 1 dip 1 on both sides of (|2.7[) . 



4 



The geometry 



The consistency relation (didjfi k — djdiHk = 0) of equation (I2.10p is given by 

T^jFik = VjFjfc + \p,j , F ik ] — Q^F mk — Q™jF im , (2-12) 

where we have used the Bianchi identities Vft-Fy + V jF ki + ViFj k = 0. Here Tl n j ik is the 
generalised curvature tensor and is defined by 

TZ" jik = dffikj ~ dk^ij + ^im^Tj ~ ^km^ij • (2-13) 

We notice immediately that if F^ = (that is, «j = M _1 <9jM for some invertible matrix 
M(ip)), then VJ 1 - ik [i n = 0. Since the matrices fa are assumed to be linearly independent, 
we have lZ n j ik = and the target space of the non-linear sigma model is, in this case, 
parallelisable. 

In the context of string theory, non-linear sigma models describe the propagation of 
strings in non-trivial backgrounds. The consistency of this propagation is equivalent to 
the vanishing of the beta functions of the non-linear sigma model [TUl [HJ H2] • At the one 
loop level, these beta functions (in the absence of the dilaton field) are characterised by 
the generalised Ricci tensor TZ n i = gi k 1Z n j ik . If the non-linear sigma model admits a zero 
curvature representation then its generalised Ricci tensor satisfies 

U n if i n = -V k F lk . (2.14) 

In particular, if the generalised Ricci tensor vanishes (namely, conformal invariance holds at 
the one loop level) then one has 

V k F lk = . (2.15) 

This last equation is a generalisation of the equations of motion of pure non-Abelian Yang- 
Mills gauge theory. 

It is also interesting to split equation (I2.10p into its symmetric and anti-symmetric parts. 
This yields 

= V i fj lj + V j fa-2T^ k 

Fij = --(Vifa-V^-H^,, , (2.16) 

The first equation is a gauged version of a matrix valued Killing equation. Indeed, if [an , + 
[dj , fa] = then this first equation is simply d^j + djfa — ^F^fa — 0. In this case the 
entries of fa are Killing vectors (isometries) of the metric . 

The linear system (12. 5p with Pi = a.i + fa could be made 'more symmetric' by writing 
®i — li ~ and ^ = ji + \fa, for some matrices 7t(y?). The linear system takes then the 
form 



* = 

* = . (2.17) 
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In terms of the matrices 7$, the symmetric and anti-symmetric parts of equation (12.81) give 

diHj + djfXi + [7j , fj,j) + [jj , fjtt] - 2r}j. fj, k = 

= . (2.18) 



1 



dcfj - d^i + [ji , 7^] - - fa , fij] + H {j fi k 

The advantage of working with the matrices 7$ is that the derivatives of fii do not appear in 
the anti-symmetric part of the last set of equations. 

Another geometric structure occurs when introducing the following change of variables 



2 9ij v 3 = (Qij + Qji) i /J 
OCi + Qu v l . 

In terms of the new variables v % and 5j, equation (12.101) takes the form 



IM 

Oti 



v l diQij + Qh ViV 1 + Qu VjV 1 = Q ik Q 



ij 



— F 

J- 1*, 



(2.19) 



(2.20) 



where Fij = didtj — djdn + [on , atj] and V jV l = dju 1 + ocj , u l . Notice that the left-hand side 
of this last equation is a gauged version of a matrix valued Lie derivative for the tensor Q^. 
With the new variables, 5j and z/, the linear system is given by 



d+ (-Qu u l + a,) dip* 
d + (Qkj v k + «i) <V 








(2.21) 



Equation (I2.20p . when = and when v % take values in a Lie algebra, is precisely the 
relation encountered in the context of Poisson-Lie duality and whose solution was given by 
Klimcfk and Severa in [T3J [TJ] . 

Furthermore, equation (I2.20p can be interpreted in the following way: Let us choose two 
currents such that 



r 
r -- 



w 



-1 



dip* 



(2.22) 



with v l a and w l a two field- dependent matrices and whose inverses are, respectively, (v _1 )° and 
(u^ 1 )". In terms of these currents, the equations of motion of the non-linear sigma model 
(12. ip are expressed as 

s l = Q ljW idi a + Quv l a dr - 



■;,><1< hQ U - v% (Q l3 w{) - w k b d k (Quv\)] I a I b = . (2.23) 



In addition, these currents satisfy the Bianchi identities (stemming from the identity dd(p l 
ddip* = in (T2J22])) 



B l = w l a dl a - v l a dl a + (v k a d k w[ - w k b d k v l a ) I a P = . 
In terms of these currents the Lax pair (12.51) takes the form 

* = 



(2.24) 



d + a iV l a r 
B + f3 jW 3 b I b 



(2.25) 
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The consistency condition (the zero curvature condition) of this Lax pair is 

Z = w{(3 3 di b - vla t Br + {v% (wlfy) - w k b d k (yiou) + v\w{ [a, , I a I b = . (2.26) 

Equation (12.20P amounts then to demanding that 

Z = i/£ l + a l B l , (2.27) 

where ai = —Qui' 1 + ol% and and fa = Quv + 5j. This means that the Lax pair yields the 
equations of motion of the non-linear sigma model up to terms which identically vanish. The 
importance of working with currents will show up in the rest of the paper. 

3 Known Solutions 

As stated above, all the quantities entering equation ( I2.10p are unknowns. In order to find 
some solutions, we proceed by fixing some of these unknowns. 

As a start, let us first check that this formalism reproduces the two well-known integrable 
non-linear sigma models, namely the principle chiral model and the Wess-Zumino-Witten 
model. These models are found by taking the following expressions for the matrices aj and 

m 



g l dig , /J>i = yg 1 d i g , (3.1) 



where g (if) is a Lie group element corresponding to some Lie algebra Q defined by the 
commutation relations [T a , T b ] = f^ b T c . The indices of the Lie algebra a , b , c , . . . have 

the same range as those of the target space of the sigma model i ,j ,k , We will use 

the fact that the gauge connection Ai = g~ x dig = ef ((p) T a satisfies the Bianchi identity 
diAj — djAi + [Ai , Aj] = 0. The inverses of the vielbiens e° are denoted E\ and satisfy 
e^E l b = 5 b and efE^ = 5j . Finally, the quantities x and y are two constant parameters which 
will provide the spectral parameter. We assume that x and y are different from zero. 
Injecting the expressions of on and fa in (I2.10p leads to 



E% = Kete)E k c f c ab , (3.2) 



where k = — ^ (x 2 — x + xy — \y 

The above Christoffel symbols Y\+ and torsion are those corresponding to the following 



metric gij and anti-symmetric tensor b. 



9ij Vab&i &j 



Hijk — K Vdafbc e i e j e k ) (3-3) 

where rj ab is an invertible bilinear form of the Lie algebra Q satisfying r\ ab j h cd + g c bfad = 0- 
Owing to the property that fte" — djt\ + fb c e i e< j = 0; the torsion is a closed three form. 
Therefore fey exists locally. 
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To summarise, the Lax pair construction for the class of theories represented by the 
metric and the torsion in (13. 3p is given by 



d + x (g 1 d i gj dip 1 
2x + 2n-l V y jy J y 








(3.4) 



with x being the spectral parameter and k a parameter defining the different models. The 
class of non-linear sigma models defined by (13. 3p includes the principal chiral sigma model 
(k — 0); the Wess-Zumino-Witten model (k = |); and the non-conformally invariant Wess- 
Zumino-Witten model (k ^ \). 

Another interesting theory is found when the matrices «j and pi are constant. In this 
case we take 



x Ti , Hi = yTi 



where [Tj , Tj] = ffj T^. Replacing these in equation (12. lOf) yields 



ij 



Pi, 



ij i 



where p = — - (x 2 + xy). 



These relations yield a non-linear sigma model defined by 



(3.5) 



(3.6) 



9ij — Vij 

2 

IJ — of'll.l.l,. 



(3.7) 



where rjij is the invertible bilinear form corresponding to the Lie algebra [Ti , I}] = /,*■ Tt 
(rjij must satisfy r/kif-j + rjuflj = in order for = pr/kif-j to be totally anti-symmetric). 
Therefore, the linear system for this non-linear sigma model (13. 7p is given by 



d + xTidtp* * 



i) + Tj Bp 1 



x + p 








(3. 



where x plays the role of the spectral parameter while p is a free parameter The quantum 
properties of this model, thought for a while to be the dual of the principal chiral sigma 
model [T5l [IB] , have been studied in [T7] . 

In this paper, we will present other non-linear sigma models which admit a Lax pair 
representation. Some of these models are new integrable two-dimensional theories. 



3 In fact the parameter p has no physical meaning as it can be simply absorbed by a rescaling of the field 
ip l . This operation leads to an overall factor in the Lagrangian. 
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4 Construction having a multiplicative spectral param- 
eter 



So far we have only given the conditions under which the equations of motion of a non-linear 
sigma model admit a zero curvature representation. However, as mentioned earlier these 
conditions (namely, equation f l2.10p ) do not guarantee the existence of a spectral parameter. 
For instance, the non-linear sigma model constructed in the context of Poisson-Lie T-duality 
[13| IT4"] admits a zero curvature representation but without a spectral parameter. Therefore, 
the class of non-linear sigma models enjoying a zero curvature representation is necessarily 
larger than the class of non-linear sigma models possessing the same property but with a 
spectral parameter. 

As the presence of a spectral parameter in the Lax pair is of crucial importance in 
extracting conserved quantities [U|2], we will now demand that our construction depends on 
such a spectral parameter. 

We assume, in this section, that this spectral parameter enters in a multiplicative manner. 
Namely, we consider the following ansatz for the linear system 



d + x Sj ((f) dip 1 



d + y(x) (3j (ip) dip J 








(4.1) 



where x is our arbitrary spectral parameter and y(x) is a function of x. The matrices and 
/3i are independent of x. 

The equivalence of the compatibility condition of the linear system and the equations of 
motion of the non-linear sigma model (J 7 = fiiS 1 ) results in the equality 

(yPi - xa^j ddip 1 + (ydify - xdj&i + xy [a { , /%] ) dtfdtf? = fi t (Bdcp 1 + tydy'dip*} . (4.2) 

Upon identification of the terms involving ddtp' 1 and d(p l d(pi on both sides of this equation 
we get 



ydi/3j — xdjdii + xy Si , (3j 



&ij fa 



The first equation of this set determines the matrix //j while the second leads to 



y 



(difij - ^ijPi) - x (djcii - fi-j-S/) +xy a i: /3 







(4.3) 



(4.4) 



The important point here is that this last equation should hold for any value of the spectral 
parameter x (we recall that 2,, fa and f^ - are independent of x). This requirement is fulfilled 
only ifl 



( 



a 



djSi 



1 1 Pj 
Si , & 



(4.5) 



4 Equation (|4.4|) is of the form yA — xB + xyC = for three ^-independent matrices A, B and C. There 
are, of course, various cases to be studied. We have analysed here only the case when the three matrices 
are proportional to each other (^4 = aC, B = bC). The other cases lead, in general, to known integrable 
non-linear sigma models. In the rest of the paper, we will assume that both a and b are different from zero. 
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together with 

bx 

y(x) = — — (4.6) 
x + a 

for two parameters a and b (which will define the different non-linear sigma models). It 
is easy to see that by taking Sj = = g _1 dig, a = (2k — l)/2 and b = (2k + l)/2 one 
recovers the construction already given in (13.41) . Similarly, by taking Sj = = T, (where 
[Tj , Tj] = Tfc) and a = b = p we arrive at the linear system in (I3.8p . 



The geometry 



It is instructive to point out the special features of the theories which admit the Lax con- 
struction fl4.ll) . We would like, in particular, to provide an interpretation for the two matrices 
Sj and fa of this Lax pair. In order to do this, we start by combining the two equations in 
(14. 5 p and taking the symmetric and anti-symmetric parts. This yields 



d i K i + d j K i -2VlK l 



diLj 



djLi + 2Hl K t 



: 





diLj + djLi - 2TI L\ = [Li , K 3 \ + \L 3 , K,\ 
diKj - djKi + 2E\- Li = [Li , Lj] - [K { , Kj 



where 



(4.7) 



Li 
Ki 



b (3i — a&i 



The first two equations of the set (14. 7ft can be cast in the form 
K l d l g ij +g lj d i K l + g il d j K l = Q 

K l d x b l3 + b^K 1 + bufyK 1 = - [di (Lj + b jt K l ) - dj ^ + b u K l 



(U 



(4.9) 



where K l = g^Kj. In terms of the tensor = gij + bij, these two equations combine to 
yield 



K l diQij + QijdiK 1 + QudjK 1 = - [di (Lj + bjiK 1 ) - dj (Li + b a K l 



(4.10) 



This last relation is precisely the condition needed for the non-linear sigma model (12. ip to 
enjoy the isometry symmetry [T8J [19] 



6<p* 



e AB K 



AB > 



(4.11) 



where K l AB are the entries of the matrix K l and e AB are constant infinitesimal parameters. 
We conclude that if a non-linear sigma model accepts a Lax pair representation of the form 
(14. ip then this model possesses an isometry symmetry (of course, for K l ^ 0). 
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The conserved currents corresponding to the isometry transformation ( 14. lip are 



J = (Ki — Li) dp 1 = —2a oiidip 1 
J = {K i + L i )8p> i = 2bp1 i dp i . 



(4.12) 



In terms of the two currents J and J, the Lax pair ( 14. ip is given by 

d - — j] V = 

2a \ 

x 

a + 



J 







2 (x + a) 

Moreover, by contacting both sides of equations (14.51) with dip 1 dpi , we obtain 



(4.13) 



0J + 



dJ 



J, J 
J, J 



2b A £ l 
—2aai£ l 



(4.14) 



where S l = ddp 1 + Q l i jdp l dp : ' = are the equations of motion of the non-linear sigma 
model. These relations can be cast into the form 



8J + 0J 



0J-8J + 



J, J 



2KiZ l 
2L i £ i 



(4.15) 



This last set of equations, suggests the study of three different cases: 
1)1^ = and U ^ : 

In this case the two currents J and J satisfy the equation dJ + d J = independently of the 
equations of motion of the non-linear sigma model. Therefore, the equation dJ + 3 J = 
is a Bianchi identity (a topological property) for the two currents. Furthermore, equations 
( 14. 7p reduce to 



diLj — djLi = 

dtLj + - 2T\. L l = 



2H tj Li 



[Li , Lj] 



(4.16) 



This set has a unique solution given by Li = 2aTj, where [Tj , Tj] = f^Tk and the integrable 
non-linear sigma model is the model defined in (13. 7p with the identification a = p. 
2) L { = and K t ^ : 

Here the two currents J and J are the conserved currents corresponding to the isometry 
symmetry 5p l = e AB K % AB of the non-linear sigma model and satisfy the Bianchi identity 



dJ-dJ + 



J, J 



0. In this case the set (14. 7p gives 

d i K j + d j K i -2T l i ,K l ^0 



2H l ij K l = 



diKj - djKi + [Ki , Kj] = 



(4.17) 
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This has a unique solution given by Ki = g^dig, for some Lie group element g(ip), and the 
integrable theory is the principal chiral non-linear sigma model (that is, equation (13. 3p with 
k = 0). The parameter a in (I4.13P is equal to —1/2. 
3) Ki ^ and U ^ : 

This case means that the integrable non-linear sigma model must possess the isometry sym- 
metry (14. lip whose conserved currents are J and J (namely, dJ+dJ = on shell). Moreover, 
the field strength d J — d J + J , J vanishes only when the equations of motion are obeyed 
and is no longer a Bianchi identity as in the previous case. 

As seen earlier, the first two cases are unique and lead to known integrable non-linear 
sigma models. Therefore, any new integrable non-linear sigma model (having (14 . 1 3 j) as a 
Lax pair) must fit in this third class. Notice that the conformal and non conformal WZW 
models are part of this third class. Indeed, by taking Ki = g _1 dig and Li = 2ng~ 1 dig (for 
k ^0) and injecting them into the set (14 .7p one gets precisely equations (I3.2p leading to the 
WZW models (13731) . 

It is therefore natural to ask whether other non-linear sigma models, possessing (14.131) as 
a Lax pair, exist. We will provide below an example of a non-linear sigma model, different 
from (13. 3 p and (13. 7p . which have such a property. 



An Example 



The non-linear sigma model we consider here is based on the SU(2) Lie algebra. Its 
action is given by 

S = J dzdz {drdr + B(r)5 ab dn a dn b + C{r)e abc n a dn b dn c ) , (4.18) 

where r and n a (a = 1,2,3) are the fields of the non-linear sigma model subject to the 
constraint 5 ab n a n b = 1. The SU{2) Lie algebra is [T a , T&] = e a b c T c with 6123 = 1. We shall 
make no distinction between upper and lower SU(2) indices. 

It is important to mention that this non-linear sigma model was first studied in [3] 
and later reexamined in j3]. The authors of [5] demanded that this model should be: firstly 
classically integrable and secondly it admits two commuting Kac-Moody algebras at the level 
of Poisson brackets (the Hamiltonian is quadratic in the Kac-Moody currents). Here we will 
require that the model is only classically integrable. Despite the fact that our requirement 
is less restrictive, the only solution we found is that of [3]. 

The equations of motion of this non-linear sigma model are 

= -2ddr + B'5 ab dn a dn b + C'e abc n a dn b dn c , 

= -2Bddn b - B' (drdn b + drdn b ) - 2Bn b 5 cd dn c dn d 



3Cdn a dn c + C'n a (drdn c - drdn c )] - 3Ce cde n b n c dn d dn e . (4.19) 



Here a prime denotes the derivative with respect to r. The equations of motion for the field 
n b are obtained by introducing a Lagrange multiplier for the constraint 5 ab n a n b = 1. 
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The ansatz for the SU(2) components of the two currents J = J a T a and J = J a T a , 
appearing in the linear system (|4.13p . is taken to have the form [3] 

ja = n *dr + (3(r)dn a + -i(r)t abc n b dn c 

ja = - n *§r-l3(r)dn a + j(r)e abc n b dn c (4.20) 

Our task now is the determination of the functions B(r), C(r), /3(r) and 7(7-) such that on 
shell (namely, when the equations of motion (I4.19P are obeyed), these two currents satisfy 

dJ a + dJ a =0 
dJ a -dJ a + e abc J b J c =0 (4.21) 

In practice, one extracts ddr and ddn b from the equations of motion (I4.19P and injects them 
into the above on-shell requirements (I4.2ip . 

We find that f)4.2ip holds provided that the following differential equations are obeyed^] 

-l + 0'-lC = Q 
-(l- 7 )-/3' + |i?' = 

27 - /? 2 - 7 2 - C = 

2/3(1-7) ~B' = (4.22) 

The first two equations of this system are solved by 

B = 67 

C = -b(r-p)+c , (4.23) 

where b and c are two arbitrary constants, with b different from zero. The solutions to the 
whole system are then divided in three cases: 
Case a) 

7 = 1 , = ±y/(l + b) (4.24) 

with b arbitrary. 
Case b) 

'^u^fr^ • ^-jr^hi- (425) 

[(r-d) -1] [(r-d) -1] 



5 We have used the fact that e a f, c e C£ j e = 8 a d$be ~ ^ae^bd and other relations coming from the differentiation 
of the constraint 8 a bii a n b = 1. In particular, the relation e bcd n a n b dn c dn d = e abc dn b dn c . 
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where d is a constant of integration. 
Case c) 



fe^-1 



1 + be^-^ 



7 



[1 _ 5 e r?a(r-d)j 2 _ 4 e »?a(r- 



/3 



b 2 e 2r,a{r-d) 



(4.26) 



where d is a constant of integration, a = |y(X+~fe) and r/ 2 = 1. 

In summary, the SU(2) integrable non-linear sigma model is given by 



5 



dzd^ {dr«9r + b-f(r)5 ab dn a dn b + [b (f3(r) - (r - d)) + (c - fed)] e abc n a «V«9n c } 



(4.27) 

It is clear from the expressions of the two functions /3(r) and 7(7*) that the constant d can 
be absorbed by the change of variable r — > r + d. Furthermore, the term involving the 
constant (c — bd) in the action (I4.27P is a total derivative and does not contribute to the 
equations of motion. Therefore, the integrable non-linear sigma model (I4.27p . up to some 
rescalings of the fields, is precisely the model of [3] . 



5 Generalisation of the principal chiral model 

We have so far studied only the case when the spectral parameter enters the Lax pair in 
a multiplicative manner. What can one say now about the situations when this is not the 
case? Unfortunately, there is not much that one can say about a general non-linear sigma 
model. The integrability is described by the Lax pair (12. 5 p and the master equation (I2.10p 
where one might find a spectral parameter hidden in the matrices on and However, a 
great deal can be learnt by exploring particular non-linear sigma models. 

We start by studying a theory which generalises the principal chiral non-linear sigma 
model. The classical integrability of this model has already been investigated in the literature 
[20| [21] . Its study here is for two purposes: Firstly, we would like to highlight the group 
theory structure behind the integrability requirement. Secondly, the results of this example 
will be of use in the rest of the paper. 

We consider the generalisation of the principal chiral non-linear sigma model as given by 
the action 



S(g) 



dzdzVt ab [g dg) [g dg 



- \b 



(5.1) 



where VL ab is a constant matrix and g 1 dg = {g 1 dg) a T a = A a T a , g 1 dg = [g 1 dg^j T a = 
A a T a and T a are some Lie algebra generators, [T a , T b ] = fa b T c , and g an element in the 
corresponding Lie group. The equations of motion are conveniently expressed as 



dA d + dA d ) + 



- (Q cd - fide) fab + (P'adfbc + ^dbfa 



A a A b = 0,(5.2) 



where use of the Bianchi identities 



B a = Q A a _ Q A a + fajbjp = Q 



(5.3) 
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has been made. 

The linear system, whose consistency condition is equivalent to the above equations of 
motion, could only be of the form 

(d + A a P a ) m = 

(8 + A b Q b )y = , (5.4) 

where P a and Q a are constant matrices. The spectral parameter, if it exists, is hidden in 
these matrices. The compatibility condition of this linear system is 

T = \ (Q a - P a ) (dA a + BA a ) + (-±/* (Q d + P d ) + [P a , Q 6 ]) A a A b = . (5.5) 

We then demand that 

7 = R a £ a (5.6) 

for some linearly independent matrices R a . This leads firstly to the determination of the 
matrices P a in terms of Q a and R a through 

Pa = Qa + (tt ac + n ca )R C . (5.7) 

Secondly, it yields the condition 

[Qa , Q b ] = f d ab Q d + si cd f d ab R c + {n ad fi + n db fi) r c - (n ac + n ca ) [r c , Q b ] (5.8) 

This last equation determines the non-linear sigma model (that is Q ab ) and its Lax pair (the 
matrices Q a and R a ). It is clear that if a solution exists, then it must have a Lie algebra 
interpretation. We will show below that equation ( 15. 81) does have solutions for particular 
examples. 



A three-dimensional example 



This particular solution to (15. 8p is based on the SU(2) Lie algebra. We take the matrix 
Q ab to be diagonal 

tt ab = diag (Li , L 2 , L 3 ) . (5.9) 
The linear system is taken to have the form 

(9 + X^A b T a ) * = 

(B + Y d A d T^ * = , (5.10) 

where the SU(2) Lie algebra is [T a , T&] = e abc T c with 6123 = 1 and and F b a are constant 

3x3 matrices. In the notation of ( 15. 4p . we have P a = X b T b , Q a = Y b T b and R d = 

l(n-y%x b a -Y a b )T b . 
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It is then found that equation (15. 8p is satisfied for the following non-vanishing elements 



X% and Y b a : 



xl 



Y 2 

1 2 



Y 3 

r 3 



1 



L2VL3 
1 

1 

/ Lia/L2 l 

1 
1 

1 



KiJ X + L 2 J X + L 3 + uiiy/xJx + Li 



k 2 \ x + Li\ x + L 3 + oj 2 \fx\ x + L 



K 3 Jx + Li\ x + L 3 + u 3 y/xJx + L 



KlA/X + L 2 \/x + L 3 - uiyfx-Jx + Li 



k 2 J x + LiJ x + L 3 - u 2 \fxJx + L 



k 3 a/x + Li Jx + L 3 - u 3 \fxJx + L 



(5.11) 



where k\ — k 2 , — 1, k 3 = kik 2 and Wi = k 2 oj 3) oj 2 = K\0J 3 , u 2 = 1. The spectral parameter is 
x. This solution was first found in [6] and later given in its present form in [2Tj . 

This solution might as well be extended to the case of the SO(6) Lie algebra which is 
described by the commutations relations 



[J a , Jb] = tabcJc , [J a , M b ] = 6 abc M c , [M a , M b ] = € abc J c 

However, this Lie algebra can be cast in the form 



°a 1 °b 



J+ , J b ~] = with J^ = -(J a ± M a 



(5.12) 



(5.13) 



We have therefore two copies of the SU(2) Lie algebras and we conclude that the non-linear 
sigma model (15. ip defined by 



n 



ab 



diag (Lf , L\ , L\ , L l , L 2 , L P 



(5.14) 



is integrable. The corresponding linear system is of the form (I5.10p where now the indices 
a , b . . . run from 1 to 6. The matrix elements (Xl , X 2 , Xf ) and (Xf , Xf , X§) are obtained 
from (15.1 ip by replacing (Li , L 2 , L 3 ) by (L^ , , L 3 ) and (Lj~ , L 2 , L3 ), respectively. 



The determination of Y b a is similar. 



A four- dimensional example 



So far, all the integrable non-linear sigma models listed in this paper are known. We 
will now explicitly construct a new integrable model. This is based on a four- dimensional 
non-semi-simple Lie algebra whose commutation relations are 

■I- I'; <,,!', , [P i ,P j \=e ij T, , [T , J] = [T , P { ] = 0, (5.15) 
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with ei2 = 1. We denote the generators of this Lie algebra as T a = {Pi , P 2 , J , T} such 
that the Lie algebra is [T a , T b ] = /„ 6 T C . Although this algebra is non- semi- simple, it has 
nevertheless a non-degenerate bilinear form r\ ab obeying rj ab f^ d + r\ cb f h ad = 0. It is given by 
[221 



Vab 



( 1 \ 

10 

6 1 

V o o i o j 



(5.16) 



where h is an arbitrary constant. 

Using this invertible bilinear form one can, for instance, construct a principal chiral 
non-linear sigma model 



S(g) 



dzdzr) ab (g dg) (g dg 



(5.17) 



where as usual (g^dg) = {g- l dg) a T a = A a T a and (g~ l dg) = (g^dgY T a = A a T a . This 
an integrable non-linear sigma model with the Lax pair 



is 



d + 



9 dg 



1 + x 

B + (g-'Bg) 








(5.18) 



where x is the spectral parameter. 

Our aim is to generalise the above principal chiral non-linear sigma model. We consider 
the action 

(5.19) 



S(g) = J dzdzQ ab (g- l dg) a (g- l dg 



where Q ab is symmetric and of the form 



n 



ab 



( Lx \ 

L 2 

b L 3 

V L 3 j 



(5.20) 



This non-linear sigma model is integrable if a solution to (15. 8p can be found. For this purpose, 
we seek a linear system of the form 



(d + X^A b T a ) * = 
(<9 + Y^A d T c ?j * = 



(5.21) 



where we assume that the constant matrices and y d c are diagonal (that is, eight un- 
knowns). It turns out that (I5.8P yields six independent equations (for eight unknowns). 
These are 

- 2L 3 Y 2 2 Xl + (Lx - L 2 ) (X 4 4 - F 4 4 ) + L 3 (X 4 4 + F 4 4 ) = 
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+2L 3 Y 1 1 X* + (L, - L 2 ) (Xt - Y*) - L 3 (X 4 4 + 
+2L 2 Y 3 3 Xl + (L 3 - LO (Xl - Y 2 ) - L 2 (x 2 2 + Y' 2 
+2L 2 Y 1 1 Jff + (L x - L 3 ) (x| - Y 2 2 ) - L 2 (x 2 2 + Y\ 



-2L 1 K S 3 X 2 2 + (L 2 - L 3 ) (X* - y/J + Li (X x x + 
-2L 1 Y 2 Xl + (L 3 - L 2 ) (x 1 - Y"/) + L a (x\ + Y"/ 











(5.22) 



The explicit solution to this set is too untidy to write down. We choose therefore to give it 
in an iterative manner. 

We take the two elements X\ and Y"/ as our free parameters leaving us with six equations 
and six unknowns. One of these parameters reflects simply the fact that our Lie algebra is 
invariant under the rescaling Pi — > XP±, P 2 — > XP 2 , J — > J, T — > X 2 T. Therefore, only one 
parameter is left and in order to simplify slightly the equations we take X 1 = 1 + x and 

= 1 — x. The parameter x will play the role of the spectral parameter. The last two 
equations of f)5.22p yield then 



1 


1 - 




1 

xl 


1 + 



x(L 2 - L 3 ) 
Lx 

x (L 2 - L 3 ) 
Lx 



(5.23) 



while the first and the second equations give 



Y 4 
1 i 



(Li -L 2 - L 3 ) (1 - x) 2 (Lx - L 2 + L 3 ) (1 + x) 2 

2(Lx-L 2 ) 2 + 2(Lx-L 2 ) 

(Lx -L 2 + L 3 ) (1 - x) 2 (Lx - L 2 - L 3 ) (1 + x) v2 

An H —7^ In 



2(Lx-L 2 ) 



2 (Lx - L 2 



(5.24) 



Replacing these back into f )5.22p allows the determination of Y 2 



Y 2 



(L 3 — Lx — Li 
(L 3 -Lx + L 2 



-X 2 + 



2(l + x)L 2 (Lx + x(L 2 -L 3 )) 1 



Lx (L 3 -Lx + L 2 ) 



X 2 



(5.25) 



Finally, X 2 is obtained by solving the equation 



L\ (Lx - L 3 ) (X 2 ) 4 - L l [(L 3 - Lx - L 2 ) ((L 3 - L 2 ) x 2 - L,) + 4L 2 (L x - L 3 ) x] (X 2 ) 2 

a " " ' 11 2 " (5.26) 



+L 2 (1 + xf (L, +x(L 2 - L 3 )f = 



This is a quadratic equation in z = (X|) 2 whose solution is 

.2 1 



2L\ (L x - L 3 



[Lx [(L 3 -Lx- L 2 ) ((Lg - L 2 ) x 2 - Lx) + 4L 2 (Lx - L 3 ) x 



± 



Lx [(L 3 - L x - L 2 ) ((L 3 - L 2 ) x 2 - L,) + 4L 2 (L x - L 3 ) x 

1/2 



- 4L 2 L 2 (L 1 -L 3 )(l + x) 2 (L 1 + x(L 2 -L 3 )) 2 } } 



(5.27) 
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Once X\ has been determined all the other elements X% and Y£ are found. 

We should mention that this solution is not valid for the two values of the spectral 
parameter x — —1 and x = L\/{L% — L 2 ) as they lead to X\ = 0. Furthermore, there are 
special cases corresponding to particular values of L\, L 2 and L 3 which are not discussed 
here. These particular cases are such that L\ — L 2 = 0, L\ — L 3 = 0, L 3 — L\ + L 2 = or 
L 3 — L\ — L 2 = 0. 



A five- dimensional example 



We consider a five- dimensional non-semi-simple Lie algebra whose generators T a 
{Pi , P 2 , P 3 , J , T} satisfy [T a , T&] = f£ b T c with the non-vanishing commutators being 



[Pi , Pj\ — Cijk Pk + V k tkij J , [T , Pi] — V k tkij Pj 

where 6123 = 1 an v k is an arbitrary constant vector. 

This Lie algebra possesses an invertible bilinear form given by [23] 



/ 1 \ 

10 

Vab= 10 

1 

\ 1 J 



(5.28) 



(5.29) 



Of course, as in the above four-dimensional example, the principal chiral model constructed 
with this bilinear form is integrable. 

The integrable non-linear sigma model we present here is described by the action 



S(g) = J dzdzQ ab (g- l dg) a (g- l Bg 



(5.30) 



where Q a b is a slight generalisation of the bilinear form (I5.29P and is of the form 



ab 



( 1 \ 

10 

10 

L 

\ L J 



(5.31) 



This non-linear sigma model is integrable and its corresponding linear system is given byH 

(d + X a b A b T a ) $ = 

(5 + Y^A d T c ) * = , (5.32) 



6 The condition of integrability (|5.8I) leads to a long list of equations which we will not write down. Here 
we will only give the result of a computer based investigation. 
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where the constant matrices X% and Y£ are 



x a b 



1 b 



( x 
0x00 
0x0 
y 










\ 1 + L(x- 1) J 

( x 
0x0 
a; 

2x 2 - y (2x 
\ 



2x 










(5.33) 



(2x-l)-L(x-l) / 



The parameters x and y are arbitrary. 



Klimcik's solutiorE 



The most important ingredient in the construction of [24J is a linear operator acting on 

the generators of a simple Lie algebra Q. This operator is required to satisfy two relations. 
Firstly 

[RX ,RY] = R([X, Y] R ) + [X,Y] for X , Y G Q , (5.34) 

where 

= , Y] + [X , RY] . (5.35) 

Secondly, R verifies the skew-symmetry condition 

< RX , Y >g + < X , RY > g = , (5.36) 

where < , >g is the Killing-Cartan form on the Lie algebra Q. 

An R operator satisfying f !5.34j) and (I5.36P was given in [21]. This is constructed as 
follows: Let the Lie algebra Q be generated by (H^ , B a , C a ) where a labels the positive 
roots and are the generators of the Cartan subalgebra. The step generators E a and E~ a 
are defined, up to a normalisation factor, as 

E a = B a + iC a , E~ a = B a -iC a . (5.37) 

The linear operator R acts on the generators of the Lie algebra Q as follows: 

RH» = , RB a = C a , RC a = -B a . (5.38) 

7 This solution was published during the reviewing process of the present paper. 
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For instance, in the case of the Lie algebra SU(2), with generators T = {7\ , T 2 , T 3 }, we 
have 



RT 



Similarly, for Q = SU (3) the operator R acts as 








( T 2 












V o 



(5.39) 



RT\ = T 2 , RT 2 = —Ti , RT 3 = , RT A = T 5 , 
RT5 = — T 4 , i?T 6 = T 7 , i?T 7 = — T 6 , i?T 8 = . 



(5.40) 



We will show below that (at least for the SU (2) case) the above action of the operator R is 
not the most general one. 

The integrable non-linear sigma model constructed in (23] is described by the action 



S (g) = J dzdz < g- x dg , (/ - eR)' 1 g~ l dg > g 



(5.41) 



where g is the group element corresponding to the Lie algebra Q, I is the identity operation 
on the generators of Q and e is a free parameter. The equations of motion of this theory can 
be cast into the form 



J, J 



R 



o , 



£ = d J - d J + e 
where the two currents J and J are defined as 

J=(I + eR)- 1 g- 1 dg , J — — (I — eR)~ x g~ x dg 
These two currents satisfy the Bianchi identity 



(5.42) 



(5.43) 



B = 



(dJ + df) + eR (dJ -dJ + e[j, J 



e 2 -l 



J, J 



RJ, J 



+ e 



J, RJ 



(5.44) 



These are found by writing g x dg = (I + eR) J and g x dg = — (I — eR) J and demanding 



that the identity d (g 1 dg^j — d (g l dg) 



Finally, the Lax pair corresponding to this non-linear sigma is given by 



9 l dg , g l dg 



holds. 



d - [e 2 - eR - 



d 



e 2 + eR 



l + e 2 
1 + x 
l + e 2 



J 



J 



x 








(5.45) 



with x being the spectral parameter. The zero curvature condition of this linear system is 

T = _ x (1 + e 2 ) (i + xeR) £ + 1 + x 2 e 2 B (5 4g) 



1 -x 2 



1-x 2 



We see that the second term vanishes identically while the first is equivalent to the equations 
of motion. 
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Of course the action (I5.4ip is a particular case of the generalised chiral principal model 
( 15. ip . For example, in the case of the Lie algebra SU(2), the model (15.411) leads to 



n 



ah 



1 







l+e 2 



(5.47) 



The anti-symmetric part of Q a b yields a total derivative in the action and one is left with the 
integrable model of Cherednik [6]. For completeness, we also give the matrices P a = X^Tb 
and Q a = Y^T b appearing in the linear system (15. 4p by comparison with (I5.45P 



1 + x 



f 1 xe 
-xe 1 
V l-xe 2 



Y 



x 



1 xe 
-xe 1 
1 + xe 2 



(5.48) 



where the lower index labels the rows. 

Similarly, the tensor Q a b corresponding to the SU(3) case is given by 



n 



ab 





I 1 


— e 
























e 


1 






























l + e 2 



















1 











1 


— e 













1 + e 2 











e 


1 






























1 


— e 
























e 


1 









Vo 




















l + e 2 


) 



(5.49) 



Here the anti-symmetric part does not result in a total derivative in the action. The explicit 
form of the SU(3) matrices P a = X^Tb and Q a = Y^Tb of the linear system (15 .4p are 



Xt 



x 



( 


1 


xe 
























— xe 


1 






























1 - xe 2 






























1 


xe 
























— xe 


1 






























1 


xe 
























— xe 


1 






























1 - xe 2 


) 



(5.50) 



The matrix Y^ is obtained from the expression of X^ by replacing x with — x and e with — e. 

As mentioned above, it seems that the action of the operator R as given in (I5.38P is not 
the most general one. Indeed, for the case of the Lie algebra SU(2) we have found that the 
most general linear operator R, satisfying the two conditions (I5.34p and (I5.36p . is given by 



RT 




(5.51) 
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(5.52) 



where a and b are two arbitrary parameters and c = a/1 — a 2 — b 2 . In this case, the action 
f)5.4ip yields an integrable non-linear sigma model of the form (15. ip with 

^ / (1 + e 2 b 2 ) —e(a + ebc) e (eab — c) 

~^ £ \£(£a6 + c) e(b — eac) (l + e 2 a 2 ) 

The corresponding Lax pair as read from (I5.45P results in a linear system of the form (15.41) 
with P a = X b a T b and Q a = Y*T b , where 

(1 — xe 2 b 2 ) xe (a + ebc) —xe (eab — c) \ 
-xe (a — ebc) 1 — xe 2 c 2 xe (b + eac) . (5.53) 
-xe (eab + c) — xe (6 — eac) (1 — xe 2 a 2 ) J 

The matrix Y£ is obtained from the expression of X h a by replacing x with — x and e with — e. 

Finally, we should also mention that the integrability of the model (15.411) holds even when 
the Lie algebra Q is non- semi- simple provided that the bilinear form < , >g is invertible (the 
invertibility of the bilinear form is used in deriving the equations of motion). This is so 
because the proof of the integrability of the non-linear sigma model (15.411) relies on the two 
relations (I5.34p and (15.361) . satisfied by the R operator, and on the existence of an invertible 
bilinear form on the Lie algebra Q. As an example, the non-semi-simple Lie algebra (I5.15P 
with the invertible bilinear form (I5.16P possesses an R operator of the fornj§ 



1 

1 + x 



RT 



( 



\ 





-1 
a 




1 

d 




-a 
-d 
-be 
-c J 



P2 
J 

V T J 



(5.54) 



where a, c and d are arbitrary parameters (b is the parameter already appearing in the 
invariant bilinear form (I5.16P ). The resulting integrable non-linear sigma model, as written 
in (15.11) . has a matrix Q ab given by 



n 



ab 



V 



1 

e 

e(a+ed) 
' (1+ec) 





e(a—ed) 
e(d+ea) 



e(ea—d) 
(1+ec) 





b+s 2 (b- 



(1-£ 2 C 2 ) 



(1-ec) 






l+£ 2 

(1+ec) 





The Lax pair of this integrable non-linear sigma model is of the form (15 .4p with P a 
and Q a = Y^Tb and where 



(5.55) 



X b a T b 



1 + x 



1 

— xe 

xe{a+ed) 
(1+ec) 





xe 
1 

xe(d-ea) 
(1+ec) 








l+ec+a:e(c— e) 
(1+ec) 





xe(ed-a) 
— xe(d+ea) 

xe[bc(l+e 2 )-e(a 2 +d 2 ) 

(l-e 2 c 2 ) 
1— ec— xe{c+e) 
(T^ecl 



(5.56) 



/ 



with the lower index counting the rows. The matrix Y" a 6 is obtained by substituting x with 
— x and £ with — e in the expression of X^. 



There are other R operators for the algebra (|5.15l) . We have chosen to write down the simplest of them. 
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6 Duality and integrability 



In this section we will put forward a method for constructing new integrable non-linear sigma 
models starting from already integrable ones. This is based on the concept of T-duality [25] 
that certain two-dimensional non-linear sigma have. 

In order to illustrate the important role of T-duality in the construction of integrable 
non-linear sigma models, we will start by considering a simple example. The discussion of 
the general case will be treated somewhere else [26J (Abelian T-duality in the context of 
integrability has been used in [5] in a very particular non-linear sigma model). We will 
consider the non-linear sigma model of the previous section whose action is 



S(g) = J &z&zn ab {g- l d g y {g- l dg) b , (6.1) 

This theory has the symmetry transformation g — > Lg, where L is a constant element of the 
Lie group to which g belongs. The duality transformations are found by first gauging this 
symmetry and adding a Lagrange multiplier term which constrains the gauge field strength 
to vanish (27J [281 EH] . This results in the first order action 

5a (g, B, B, X ) = J dzdz \n ab (g- l dg + g^Bgf (g^Bg + g^Bgf 

+ X a (dB a - BB a + f b a c B b B c ) ] . (6.2) 



Here B = B a T a and B = B a T a are the gauge fields with the gauge transformation B — > 
LBL^ 1 — dLL^ 1 , B — > LBL^ 1 — dLL^ 1 . The Lagrange multiplier Xa, transforming in the 
adjoint representation, imposes the pure gauge condition dB a — dB a + f£ c B b B c = whose 
solution is B = h~ 1 dh and B = hr x dh. Upon replacing this back into the action S\ one 
gets the relation S\ = S(hg) and by choosing h = 1 (thanks to the local gauge symmetry 
g — > Lg, h — > hL~ v ) one concludes that the action S\ is equivalent to the action S. 

The dual action is obtained by keeping the Lagrange multiplier and eliminating, instead, 
the gauge fields (through their equations of motion). This procedure yields, after the gauge 
choice g — 1, the dual theory 

Six) = I ^dz (M- 1 ) ab dxadxb 

M ab = n ab + Xc f c ab ■ (6.3) 

We will show now that if the original theory ( 16. ip is integrable then its dual (I6.3P is also 
integrable. It is convenient, for this purpose, to introduce the two currents 

J a = (M- l ) ba d Xb , J a = - (M- l ) ab B Xb . (6.4) 

In terms of these, the equations of motion of the dual theory ( 16. 3 j) are 

£ a = d.J a - d.J a + ft c J b J c = . (6.5) 

We notice that these are the Bianchi identities ( 15. 3 p of the original theory with ( J a , J a ) inter- 
changed with (A a , A a ). Furthermore, these currents satisfy the Bianchi identity (stemming 
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from ddxa 

B r = 



1 

~2 



= 0) 



8J d + dJ d ) + 



(Med - &<dc) fab + [p-adfbc + ^dbf, 



J a J b 



(dJ d - Bj d + fiJ b r) = o 



(6.6) 



Again these are a linear combination of the equations of motion ( 15.21) and the Bianchi 
identities (15.3p of the original sigma model with ( J a , J a ) interchanged with (A a , A a ). 
Therefore, for the dual non-linear sigma model, the linear combination 



B c + 



2 i^cd 



2 (^cd + ^dc) 



V d c) + XaF cd 

dJ d + 3 J d 



£ 



+ 



- (&>cd — &dc) fab + {P'adfbc + &>dbf t 



J a J b (6.7) 



takes exactly the form of the equations of motion (15. 21) of the original theory with the 
exchange ( J a , J a ) i — > (A a , A a ). 

As already shown in equation ( 12.271) . the zero curvature condition stemming from a Lax 
pair is a linear combination of the equations of motion of the sigma model and some corre- 
sponding Bianchi identities. Moreover, the equations of motion and the Bianchi identities of 
the dual theory are simply linear combinations of those of the original theory with the ex- 
change (J a , J a ) < — > (A a , A a ). Consequently, if the original non-linear sigma model (16. ip is 
integrable (namely, if equation (15. 8p is satisfied) then its dual theory (I6.3P is also integrable. 
The Lax pair of the dual theory is 

(d + J a P a ) ^ = 

d + J b Qb)v = . (6.8) 

This is simply the Lax pair of the original theory with ( J a , J a ) and (A a , A a ) interchanged. 
Indeed, if the original non-linear sigma model ( 16. ip is integrable (that is when (I5.8P is satis- 
fied) then the zero curvature condition of this last linear system is 



Z = R a B a + Qa + (n ba + Xcfba) ^ 



(6.9) 



where we have used the relation P a = Q a + (fi ac + Q C a)R c and the integrability condition 
(15. 8p . We notice that the first term vanishes identically and the second term yields the 
equations of motion of the dual non-linear sigma model. 



Examples 

1) The first example which enters into the class of non-linear sigma models (16.11) is of 
course the principal chiral non-linear sigma model for which the matrix Q a b — Vab, where r] a b 
is the invariant bilinear form of the underlying Lie algebra (that is, rj a bf b d + Vcbfad = 0)- fts 
corresponding dual non-linear sigma model is written in (I6.3P with 

M ab = Vab + Xcf C ab ■ (6-10) 
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The Lax pair of the dual of the principal chiral non-linear sigma model is read from (13. 4p . 
for k = 0, according to the above prescription. This is therefore given by 



d + x 



8 + 



x 



2x 



(M- 1 )" d X a 

(m- 1 )"" 8 Xd 








(6.11) 



where x is the spectral parameter. 

2) The second example we consider is the SU(2)-based non-linear sigma model whose 
corresponding matrix Q ab is given in (15.91) and for which 

M ab = S ab L b + Xc tcab • (6.12) 
Its dual partner is given by the action (16. 3p with the matrix M _1 explicitly given by 



( xl + UU 



D 



X1X2 - X3L3 X1X3 + X2L2 



X1X2 + X3L3 x{ + LiL 3 



X2X3 ~ X1L1 



(6.13) 



V X1X3 - X2L2 X2X3 + X1L1 xl + LiL 2 
where D = L1L2L3 + L\x\ + -^2X2 + ^3x1- The Lax pair of this dual theory is then given by 

[M- 1 ) eb d X e}T a }v = 

9Xf 



d + X* 



8 + Y d c 



1 \ df 



T c * = 



(6.14) 



where T a are the generators of the SU(2) Lie algebra and the matrices and Y b a are as 
given in (15.111) . 



7 Conclusion 

The question of classical integrability of two-dimensional non-linear sigma models has been 
addressed in this paper. We have first focused on the issue of representing the equations 
of motion of the non-linear sigma model as a zero curvature condition of a linear system. 
This is regardless of whether the Lax pair depended or not on a spectral parameter. This 
requirement resulted in a master equation with some interesting geometrical properties. In 
particular, it is shown that in the case when the matrices involved in this equation are Lie 
algebra valued matrices, this master equation is a generalisation of an equation encountered 
in the context of Poisson-Lie T-duality. It is therefore hopeful that a general solution along 
the lines of [T3l [H] might be found to this master equation. 

We have then put special emphasis on those constructions admitting a spectral parameter. 
Two situations emerged from this analyses. The first consists of those constructions where the 
spectral parameter enters in a multiplicative way. The geometry of the integrable non-linear 
sigma models is in this case tractable. The isometry symmetry of these sigma model plays 
an essential role and is responsible for their integrability. Even in this simplified version, the 
general solution to the master equation remains a challenging problem. The second situation 
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concerns those constructions for which the spectral parameter does not manifest itself in a 
multiplicative manner. Here, we have studied only the integrability of a generalisation of the 
principal chiral non-linear sigma model. The master equation for these models, (j5.8p . has a 
certain Lie algebra structure and one hopes that this might be helpful in finding solutions. 
We have carried out a computer assisted study and found two new integrable non-linear 
sigma models. These are given in (I5.20p and (15.311) and are based on two non-semi-simple 
Lie algebras of dimension four and five, respectively. 

Always in the context of the generalised principal chiral non-linear sigma model, we have 
given a brief summary of a new integrable non-linear sigma model that has recently been 
found by Klimcfk [24] . This solution to the master equation holds for any simple Lie algebra 
and relies on an R operator which acts on the generators of the Lie algebra. We have pointed 
out that the action of the R operator as given in [21], equation H5 .381) . is not the most general 
one. We have worked out the most general R operator for the case of the SU(2) Lie algebra. 
This has led to a more general integrable non-linear sigma model (I5.52p in comparison with 
the one in (I5.47p . However, the problem of constructing the most general R operator for 
other Lie algebras is still an open issue. Furthermore, it is shown that the R operator can 
also be extended to non- semi- simple Lie algebras. This was carried out for a particular 
example resulting in another integrable non- linear sigma model (15. 55 p . 

We have also shown that there is a connection between T-duality and integrability of 
non-linear sigma models. More precisely, if a non-linear sigma model is integrable and 
admits a T-duality transformation then its dual is also integrable. This might not sound 
surprising as two non-linear sigma models related by a T-duality transformation are by 
definition equivalent. However, it is not obvious how to find the Lax pair associated to the 
dual theory starting from the Lax pair of the original non-linear sigma model. We have 
given here the recipe for this passage. T-duality is therefore a mean for constructing new 
integrable non-linear sigma models. 

Among the open problems that could be addressed in the light of this work would be the 
extension of the standard procedure of the dressing transformations [3TJI [Tj, [2] encountered 
in the principal chiral sigma model [3"T| I3"2"j 155] to the new integrable models found in this 
paper. The construction of the conserved charges is also another challenging issue. 

Finally, we should mention that the study of the integrability of non-linear sigma mod- 
els carried out here could be of interest to string theory in its quest for integrable string 
backgrounds [3U E2J EE] • 
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